Does U — S — ex Form a Quasi-Bound State? 

Htun Htun Oo^, Khin Swe Myint^, Hiroyuki Kamada^ and Walter Glockle^ 

"^ Department of Physics, Mandalay University, Mandalay, Myanmar 
^ Department of Physics, Faculty of Engineering, Kyushyu Institute of Technology, 

Kitakyushyu 804-8550, Japan 
^^ ■ ^ Institute for theoretical Physics II, Ruhr- University Bochum, D-44 "^80 Bochum, 

f^ , Germany 

(N 

pH , We have investigated the possible existence of a quasi-bound state for the S — S — a 

r-{ ' system in the framework of Faddeev calculations. We are particularly interested in the state 

of total iso-spin T=2, since for an inert a particle there is no strong conversion to S ~ N ~ a 
or yl — yl — Q possible. A E — a optical potential based on Nijmegen model D and original 
S — S interactions of the series of Nijmegen potentials NSC97 as well a simulated Gaussian 
type versions thereof are used. Our investigation of the S — E — a system leads to a quasi 
bound state where, depending on the potential parameters, the energy ranges between -1.4 
and -2.4 MeV and the level width is about 0.2MeV . 
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^^ • §1. Introduction 

f^ , Strangeness S=-2 hypernuclei provide information on baryon-baryon forces in 

^T ■ tiie state of S=-2. Only tliree nuclei have been identified so far, ]^j^Be, ^y^H^ and 

^^B^ The challenge is to understand their binding energies and decay properties. 

These nuclei are especially interesting since the S=-2 two-baryon system is rich in 

structure due to the conversions between AA , SN and EU. Baryon-baryon forces 

^ . for S=0,-1 and -2 are being investigated in the meson exchange picture5*''5l'Ell2J or 

r^ [ using quark models.'^ While there is a wealth of data for S=0, which allows to 

fix force parameters, the situation is still much open in the S=-l and -2 sectors. 

Therefore additional information is needed. In this study we would like to focus on 

r> I the system S — S — a in the state of total iso-spin T=2. If the a-particle would be 

C^ ■ inert, that system could not convert ioS — N — aorA — A — a. Therefore in case 

the forces would be strong enough, there might exist a low lying state with a small 

width. The width would be caused by yl — Z" conversion leading , for instance, to 

Z! + A + N-\- {^He/^H), where {N +^ He/^H) is in a state of total iso-spin T=l or 

into U + A+'^He* (T=l). 

We investigate that system S — U — a under effective simplifying assumptions. 
The U — a interaction is modeled via an optical potential based on the Nijmegen 
model D and the E — U interaction in the state of total iso-spin T=2 is taken either 
directly as a meson theoretical Nijmegen potential of the type NSCO'ii'^ or a simulated 
version thereof of the Gaussian type.^ That 3-body system is solved precisely in the 
Faddeev scheme. 

In section |21 we derive the Faddeev equations for the X! — U — a system. The 
numerical results are displayed in section |31 The Appendix shows technical details 
of the formulation. We end with a brief summary of our work and an outlook in 
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section 0] 

§2. The Faddeev Equations for the System H — U — ex 

We assume the existence of a quasi-bound state where the width of that state is 
caused by the (absorptive) imaginary part of the effective U — a potential. If the a 
particle would be inert and the two i7's couple to iso-spin T=2, the system cannot 
convert into S — N — aoiA — A — a. We focus on such a state with T=2. The 
Schrodinger equation reads 

{Ho + VeE + VE,a + Ve.o) ^ = E^ (2-1) 

It is convenient to number the two i7's as 1 and 2 and the a particle as 3. Then 
the Schrodinger equation converted into an integral equation reads 

•^ = TT^TT (^12 + ^13 + V2Z) ^ (2-2) 

-C/ — -no 

which suggests the decomposition 

1^ = -012 + V-IS + ^^23 (2-3) 

with 

Because of the identity of the Z"s and the antisymmetry of ^ the two Faddeev 
components ^23 and '(/'13 are not independent and are related to each other by 

^^13 = -Pl2^2Z (2-5) 

This then leads to two coupled equations 

V'12 = ^ ^ ^ ^12 {ipi2 + (1 - P12) V'23) (2-6) 

tj — Hq 

■023 = -^ 7r^23 (--Pl2V'23 + V'12) (2-7) 

-c/ — -no 

In a standard manner one introduces the two-body t-matrices 

tl2 = V^12 + 1^2 ^ Tr^l2 (2-8) 

-C/ — -flQ 

t23 = "^^23 + "^^23-^ Fr*23 (2-9) 

-C/ — -Hq 

and obtains the final set of two-coupled Faddeev equations 

V'12 = ^ Tr^l2 (1 - P12) 023 (2-10) 
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1p23 = ^ fr*23 {-Pl2i^23 + ^u) (2-11) 

-C/ — -no 

The total state is then given as 

<Z^=(l-Pi2)V23 + ^12 (2-12) 

We solve that system in momentum space and use a partial wave decomposition. 
To that aim we introduce two types of Jacobi momenta in terms of the individual 
momenta kj , i=l,2,3; 

P3 = ^(ki-k2) (2-13) 

q3 = 2mi:k3-m^(ki + k2) ^^.M) 

and 

mak2-ms'ki 
Pi = ] 2-15 

(mo, + mi;) ki - mi: (k2 + ks) /„ -,„n 

qi = -T. (,2-ibj 

rUa + 277117 

The momenta pa ,q3 and pi ,qi are adequate for the Faddeev amplitudes ■012 and 
-023 which are driven by the two-body t-matrices ti2 and t23, respectively. Related 
to these momenta are partial wave projected basis states 

b3?3 {hsi) J3A3 O3A3) JM(11)2) = \p-ma:i) (2-17) 

and 

Piqi 64) h (ai0 h ijih) JM(11)2\ = \piqiai) (2-18) 

The sequence of discrete quantum numbers denote orbital and spin angular 
momenta, their intermediate couplings and the couplings to the total 3-body angular 
momentum J with magnetic quantum number M. Finally (11)2 denotes the iso-spin 
coupling. Because of the identity of particles 1 and 2, {I3 + S3) has to be even. This 
imposes the only restriction to these intermediate quantum numbers. 

By a standard procedureP^t the coupled set of equations (|2-lUI2-lljl is projected 
onto those basis states. It results 
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i'u ipsqs) = {psqsas 1012) = -2 -2- Y] / dqiqf 

^ 2/^3 2M3 "1 

' f ql 

dxti2 \p^,'K2,E- -f- 

1 V 2M3, 

Ga3aifegi^)V'2lK9l) (2-19) 
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^ 2^1 2Mi "3 '^ 
-^ 2^(1 2Mi a'l 

Ga,a; Mxx) ^il i<q'i) (2-20) 

The purely geometrical quantities Gaa' resulting from recouplings are given in 
the Appendix. Further the shifted arguments in the two-body t-matrices and the 
Faddeev amplitudes under the integrals are given as 

vr2 = \ /gi + pjqj + 2piqiqsx tts = y/ piqj + q'j + 2p2qiqzx 

■k'^ = \ Jpiq{ + gj + 2p2gig3a : vr^ = ^ q{ + p{qi^ + Ip^qiq^x 

71-4 = \l p^qi + gf + 2pqiq[x vr^ = a/^i + p^q'i + 2pqiq[x 

where pi = \,p2 = w^+M^ and p = m^+m^ 

The reduced masses are 

P3 = \Me 

,^ _ 2McMs 
^•^3 — 2Ms+Ma 



Pi 



Ma+Ms 



,^ _ Ms{Mjj+Ma) 

This is an infinite set of homogenous coupled integral equations in two variables, 
which can be truncated since the two-body t-matrices drop quickly in magnitude with 
increasing angular momenta. 

§3. Results 

The set of coupled equations (|2-19II2-2U|1 are discretized in a standard manner. 
We choose Gaussian quadrature points in the variables q and x and spline interpola- 
tion for the variables vr under the integrals. The two-body t-matrices are generated 
from the Lippmann Schwinger equation again using Gaussian quadrature discretiza- 
tion. We refer for numerical details to.l^'l^ The energy eigenvalue E is determined 
as follows. The homogenous set of coupled equations is schematically written as 

r]{E)i> = K{E)i} 

with r]{E) = 1 at the energy eigenvalue. Without knowing E, one first determines 
r] and varies E such that finally rj{E) = 1. The eigenvalue rj is determined either by 
a simple power method or by a Lanczos type algorithm. For details see.1I''1^ 

In order to demonstrate our numerical accuracy we would firstly like to display 
results on the AAa system, where for model forces we recalculated ^aHc binding 
energies. The model forces are of Gaussian types and we refer to.l^*^ That sys- 
tem has been investigated before by Filikhin et al)^ using Faddeev equations in 
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Iac 


Iaa 


Ours 


Filikhin 


Myint 








-6.880 


-6.880 


-6.880 


0,1 





-6.986 


-6.987 


-6.983 


0,1,2 


0,2 


-7.050 


-7.045 


-7.041 


0,1,2,3 


0,2 


-7.084 


-7.078 


-7.073 



Table I. Comparison of \y[He binding energies in MeV based on the model forces given ills' ancP^ 
for an increasing number of partial waves. 





Vi (MeV) 


V2 (MeV) 


NSC97a "sim" 


5274.2576 


-292.7193 


NSC97C "sim" 


4886.7830 


-289.8740 


NSC97e "sim" 


4588.5040 


-294.2548 



Table II. The potential strength parameters of the Gaussian form "sim" given in Ea. l3TII for the 
S — S system. The range parameters are fii = 0.37 fm and /i2 = 1.0 fm. 

configuration space and by Myint et aLl^ using a variational method based on a 
Gaussian expansion. Our method is mathematicahy different from theirs. We solved 
the Faddeev equations in momentum space which allows us to treat any type of two 
body forces. We show in Table ^^^iJe binding energies for an increasing number of 
relative orbital angular momenta within the Aa and AA sub-systems. Our results are 
in very good agreement with the other two methods which underlines the reliability 
of our treatment. 

Now we turn to the central topic, the UUa system in the state of total iso- 
spin T=2. For the E — E potential we use either the original Nijmegen potentials 
NSC97a,c,e^ or the simulated Gaussian forms thereof .^^' The latters are given as 



Vir) = Y,V^e 



-(r-//^^)" 



(3-1) 



with the parameters shown in Table ITll The U — a potential is chosen to be com- 
plex to provide for absorptive processes, like the ones mentioned in the introduction. 
We use the form 



V, 



Ea\ 



Y^y^e 



-ir/,.f + Y,U^e- 



ir/l.,? 



(3-2) 



with the parameters Vi = —21.3 MeV, V2 = 4.8 MeV for the real part and 
Ui = 4.07 MeV, U2 = —11.73 MeV for the imaginary part. Further, one has ^i = 1.3 
fm and ^2 = 1-7 fm. 

The S — a potential has been constructed in the following manner. An effective 
U — N potential was firstly derived from the original Nijmegen model D interaction 
in the Brueckner framework. Then that effective potential was expanded into a five- 
range Gausian form. That potential was used in a generalized Hartree-Fock method 
to generate the effective X! — a potential.-^ The imaginary part arises due to AN 
to UN conversion. 
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Fig. 1. Energy trajectory in the complex energy plane for the S — a potential. The numbers at 
each point indicate the multiplicative factors by which the attractive real part and the over all 
imaginary parts were multiplied. 



Before we present the result for the three-body system EEa we firstly investigate 
the properties of the underlying two-body sub-systems. The complex U — a potential 
leads to a complex energy eigenvalue. We locate the one with the "lowest" energy 
in the following manner. We neglect the imaginary part and multiply the attractive 
real part by some enhancement factor. Choosing for instance that factor to be 2.5 
we find a binding energy of -1.0 MeV. Next we allow the imaginary part to increase 
from in steps of 0.1 until we reach the physical value 1. In this manner we find 
the energy trajectory in the complex energy plane shown in Fig.l. We end up with 
the complex energy position in the lower energy half plane just below the unitarity 
cut from to infinity. To reach that final position we have chosen that detour in the 
complex energy plane which appears to us more feasible than a direct energy search 
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Gaussian 


Nijmegen 


NSC97a 


-2.253 


-2.250 


NSC97C 


-2.460 


-2.437 


NSC97e 


-3.214 


-3.122 



Table III. Comparison of the two-body EE binding energies in MeV for the original Nijmegen 
potentials NSC97a,c,e and the simulated Gaussian potentials. 





NSC97a 


NSC97C 


NSC97e 


Nijmegen 


-1.840 


-2.378 


-2.728 


Gaussian 


-1.841 


-2.400 


-2.819 



Table IV. The energy eigen values in MeV for the SSa system using the simulated Gaussian and 
the original Nijmegen E — S potentials together with the real part alone of the E ~ a potential. 





NSC97a 


NSC97C 


NSC97e 


Nijmegen 


-1.418- i0.202 


-2.34- i0.014 


-2.376- i0.191 


Gaussian 


-1.492- i0.218 


-2.323- i0.017 


-2.354- i0.211 



Table V. The complex energy eigen values in MeV for the EEa system using the simulated 
Gaussian and the original Nijmegen E—E potentials together with the complex E — a potentials. 



for the physical potential strength. Thus we find that the effective E — a potential 
is not strong enough to generate a complex energy with a negative real part. The 
energy search in the complex energy plane was greatly simplified by using a method 
of analytical continuation in the form of the point method.—^ 

The chosen E — U potentials support a bound state. The corresponding bind- 
ing energies are displayed in Table ITTTI for the original Nijmegen potentials and the 
simulated ones. The phase shifts in the state ^Sq and T=2 for the original and sim- 
ulated potentials agree perfectly well with each other as shown in Fig. 2. Though it 
is presumably unrealistic that two i7's are bound it is the result of the meson based 
Nijmegen potentials, whose parameters have been fixed to very many data in the 
nucleon-nucleon and hyperon-nucleon sectors and where SU(3) symmetry arguments 
allow for a prediction to the S=-2 sector. We shall comment below on an ad hoc 
weakening of those potentials. 

We are interested to see the outcome of those dynamical assumptions for the 
EEa system. In this exploratory investigation we restrict all orbital angular mo- 
menta to be zero. In the first step we neglect the imaginary part of the E — a 
potential. It turns out that in all cases using the original Nijmegen E — E potentials 
or the simulated ones, the three-body system is bound. The results are given in 
Table IIVI Then we switch on the imaginary part of the E — a potential in steps of 
0.1. The results are displayed in Fig. 3 and Fig. 4 for the two cases. The effect of the 
imaginary part in the E — a potential shifts the real part of the energy slightly to 
the right and introduces a small negative imaginary part. The resulting final energy 
positions are displayed in Table O 

In order to explore the outcome using E — E potentials which are weaker than 
the ones we used, we multiplied them by overall factors 0.9 and 0.8. The resulting 
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Fig. 2. Comparison of EE-phase shifts for the original Nijmegen potentials NSC97a,c,e and the 
Gaussian potentials for ^5*0 T = 2. 
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Fig. 3. Energy trajectories in the complex energy plane for the %^He system using the original 
E — E Nijmegen NSC97e interaction. The numbers at each point indicate the multiplicative 
factors by which the attractive real part and the over all imaginary parts were multiplied. 



X! — U binding energies for the original NSC97e potential are -1.5 and -0.4 MeV, 
respectively ( in comparison to the original value of -3.1 MeV). Then we performed 
again the energy search for the E — E — a system starting with zero imaginary part 
for the effective X! — a potential. In case of the reduction factor 0.9 (0.8) this leads 
to a 3-body binding energy of -1.08 MeV (-0.32 MeV). For the full imaginary part 
we end up with (0.68 - i 0.21) MeV for the factor 0.9 and with(-0.279 - i 2. Ox 10"^ 
) MeV for the factor 0.8. We conjecture that a low lying quasi bound state or a 
resonance for the U — U — a system in the T=2 might exist in reality. 

Finally we would like to mention a possible means to approach such a state 
experimentally. One could think of the two step process of (i^~, K^) reaction on a 
^Li target; 
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Fig. 4. Energy trajectories in the complex energy plane for the %^He system using the simulated 
Gaussian form for the E — E Nijmegen NSC97e interaction. The numbers at each point indicate 
the multiplicative factors by which the attractive real part and the overall imaginary parts were 
multiplied. 



K- +6 Li ^ ^0 + TfO + n +1 He 
TT^ + n-> S- +K+ 

In this manner one populates two Z"s together with a. It remains of course the 
task to estimate the reaction rates. 



§4. Summary and Outlook 



We developed a Faddeev code for the effective three-body systems AAa and 
X!Sa. It is formulated in momentum space and is applicable for any type of two- 
baryon forces. This allowed us to use directly the original Nijmegen forces, which is 
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the first time, to the best of our knowledge, that they have been applied for these 
systems. We tested our code in a bench mark model study for the AAa system 
reproducing perfectly well the results from two earlier studies. Our results for the 
UEa system lead to a quasi bound state with a small negative imaginary part. The 
negative real part of the energy ranges between -1.4 and -2.4 MeV. These numbers 
are based on the original or simulated Nijmegen potentials for the US system in the 
state T=2, which support a bound state with a binding energy of about -2.5 MeV. 
Further we use a X! — a optical model potential, which by itself supports a complex 
energy eigen value of about (-2 -i 0.1) MeV. We also artificially reduced the overall 
strength of the E — U potential by factors 0.9 and 0.8, which moved the 3-body 
U — U — a energies towards -0.68 and -0.279 MeV, respectively, with small widths. 

Both dynamical assumptions on the U — U and E — a potentials should be 
critically reinvestigated in the future. Upcoming meson based U — E potentials 
without a bound state should be used and in addition the effective E — a potential 
should be generated more consistently using realistic a particle wave functions in 
conjunction with i7-nucleon forces related to the same theoretical model as for the 
E — E interaction. A low lying state for the E — E — a system with isospin T=2 
would provide interesting additional information on the dynamics in the strangeness 
S=-2 sector. 
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Appendix 

The purely geometrical coefficients occuring in Eas ()2-19|) and H2-20() are 



{I' - f 

X' I /' 

L s J 

li+h=l l[+l'2=l' 



{2l + i)\ (21' + 1)1 ^f k h l\r(i\f-m^ 
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/' ^ J k h f" 



j, \\ f. \ C {f'%f;00)C {/"hf -00) 



with 



and 



gh= I dxPhix)- — -f—- T 

-1 |p2qi + q3l Iqi + piqsl 



Gaia[ {qi,qi,x) 



lEH- 



LS 



E E /■ 

ll+l2=l l[+l'2=l' 



I \ 3 
jIj'I'lXl'X' LS { X I I } { 

L S J 



[+l2 



[1l) 



'{21 + 1)1 (2/' + l 



2Zi!2Z2! V 2Z1!2/^! 



'I97'J Z^l A 



^2 / 

X L f 



A' I /' 
L S J 



C(/2A/;00) 



EC(W";00)'| J, \]^ ^^, }c(r/^/; 00) C (/%/'; 00) 



with 



ff/i 



1 



1 



dxPh{x)- 

-1 |pqi + qil |qi + pqil 



1 

2 
3 

4 

5 

6 
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